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SUMMARY 

An  attempt  is  made  to  evaluate  the  unsteady  response  of  a flat  plate  airfoil  to 
large  nonpotential  flow  disturbances  in  the  form  of  a translating  rectangular  grid  of 
eddy-array.  A suitable  streamfunction  to  represent  the  translating  nonpotential  vor- 
tex-array is  chosen.  The  method  of  analysis  adopted  is  the  singularity  distribution 
principle  in  combination  with  the  time-marching  technique.  The  problem  is  solved  in 
two  stages,  namely,  (1)  auxiliary  solution  and  (2)  time-marching  solution.  By  auxiliary 
solution  is  meant  the  solution  of  the  problem  which  completely  neglects  the  presence  of 
the  wake  vortex  sheet  and  treats  time  as  a parameter;  this  results  in  a steady  flow  type 
of  analysis.  The  time-marching  part  of  the  analysis  increments  time  by  equal  steps 
starting  from  zero  time,  makes  use  of  the  auxiliary  solution,  keeps  track  of  the  shed- 
ding and  growth  of  the  wake  vortex  sheet,  evaluates  the  unsteady  response,  and  continues 
along  the  time-axis  up  to  any  specified  maximum  time  limit.  Preliminary  numerical  re- 
sults from  a computer  program  are  presented. 


SYMBOLS 


ON  STATEMENT^ 
for  public  release; 
uion  Unlimited 


Semichord  length  f 

Lift  coefficient 
Moment  coefficient 

Eddy  intensity  parameter  

Strength  of  1th  youngest  wake  vortex 

Kernel  function  defined  in  text  copy  m 

Sides  of  rectangular  eddy  cell  e- 

Bessel  function  of  kth  order  j^|  S f 

Kernel  function  defined  in  text 
Lift  force  (+ve  up) 

Moment  about  leading  edge  (+ve  clockwise) 

Time  step  number;  summation  index 

Maximum  number  of  wake  vortices  handled  by  program 
Pressure 

Polar  coordinates 

Arc  length  along  wake 

Location  vector  associated  with  S 

Time 

Period 

x-  and  y-components  of  velocity 

Mean  translational  velocity  components  of  freestream 
Total  velocity  vector 
Complex  velocity  potential 

Mean  freestream  velocity;  wake  vortex  sheet 
Space  fixed  coordinates 

Initial  location  coordinates  for  ( £, n) -coordinates 
Complex  variable  x+iy 
Z6  Variables  defined  in  text 

Direction  of  mean  freestream  velocity 
Bound  vortex  strength 
Wake  vortex  distribution 
Bound  circulation  strength 
Strength  of  a general  vortex 
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A small  value  ( < < 1 ) 

Complex  variable  %+ir 
Inclination  of  eddy-array  with  foil 

Translating  natural  coordinates;  transformed  plane  coordinates 

Fluid  density 

Variable  defined  in  text 

Stream  function 


SUBSCRIPTS: 

A Auxiliary 

E Due  to  eddies 

1 Due  to  1th  youngest  wake  vortex 

M Mixed 

n nth  time  step 

r Radial  component 

S Due  to  steady  translation 

U Unsteady 

0 Tangential  component 

SUPERSCRIPTS: 

' Nondimensional  quantity;  dummy  variable 

1 Due  to  1th  youngest  wake  vortex 

Complex  conjugate 


BACKGROUND 

The  unsteady  aerodynamic  response  of  airfoils  due  to  a time-variant  incoming  flow 
field  leads  to  many  undesirable  results.  Some  of  these  deleterious  effects  are:  (1)  loss 
of  aerodynamic  performance  (2)  structural  failure  due  to  aerodynamic  overloading  (3)  vi- 
bration which  leads  to  other  problems  like  fatigue  and  wear  (4)  noise  due  to  the  acoustic 
dipoles  which  always  accompany  such  fluctuating  forces  (5)  laminar-turbulent  boundary 
layer  transition  (6)  cavitation  and  the  attendant  pitting  and  erosion  of  working  sur- 
faces. A second  class  of  unwanted  phenomena  originates  when  the  unsteady  response  and 
tne  resulting  direct  effects  such  as  physical  movement  of  the  foils  start  interacting 
strongly  witn  tne  incoming  flow  thus  establisning  a closed  feedback  system.  These  phe- 
nomena could  be  collectively  termed  as  loss  in  stability.  Classical  flutter,  stall 
flutter,  surge  and  rotating  stall  are  some  of  the  v/ell  known  members  of  this  group. 

Tne  type  of  time-variant  inlet  flow  dealt  with  will  depend  on  the  engineering  ap- 
plication for  tne  theory.  In  an  axial  turbomachine  the  flow  disturbances  exciting  tne 
Dlades  of  a particular  row  are  either  due  to  gross  inlet  distortion,  to  wakes  from  the 
blades  of  the  preceding  row  and  structural  members  like  struts,  or  to  turbulence.  All 
these  three  types  of  flow  disturbances  could  be  assumed  to  consist  of  eddies  having 
definite  wavelengtns.  Tne  wavelengths  associated  with  these  eddies  in  terms  of  the 
blade  chord  length  may  be  taken  to  be  very  large  for  gross  inlet  distortion,  of  unit 
order  for  wakes  from  preceding  blade  rows,  and  very  small  for  freestream  turbulence. 

Since  the  effects  of  viscosity  and  turbulence  in  fluid  flows  is  to  make  the  flow  ro- 
tational these  eddies  in  all  likelihood  are  nonpotential.  The  flow  velocities  associ- 
ated witn  them  could  range  from  very  small  to  the  order  of  the  mean  flow  velocity.  Since 
one  is  mainly  interested  in  short  durations  of  these  flows  such  as  the  time  to  travel  a 
few  chord  lengths  distance  these  eddies  may  be  taken  to  be  'frozen',  stable,  or  having 
no  mutual  interference.  In  the  present  work  an  attempt  is  made  to  formulate  a reason- 
able model  for  such  time-variant,  nonpotential,  incoming  flows  and  the  attendant  unsteady 
response  of  a single  tnin  airfoil  is  evaluated. 

PROBLEM  STATEMENT 

Consider  (Fig.  1)  the  two-dimensional,  unsteady,  rotational  flow  field  of  a 'frozen' 
rectangular  grid  of  eddies  with  alternating  signs  of  equal  circulation  strength  trans- 
lating at  a constant  velocity  U,V  with  respect  to  a system  of  space-fixed  coordinates 
x,y.  The  fluid  involved  is  assumed  to  be  incompressible  and  inviscid.  The  field  of 
velocity  components  associated  with  the  eddies  alone,  neglecting  the  constant  transla- 
tional motion,  is  represented  by  the  streamfunction 

*(?,")  « g cos~jr  C0ST 

where  i,n  is  a system  of  natural  coordinates  translating  with  the  grid  at  the  constant 
velocity  U,V;  g is  an  intensity  parameter  for  the  strength  of  the  nonpotential  eddies; 

K and  n are  parameters  representing  the  unequal  lengths  of  the  sides  of  the  rectangle 


wnicn  defines  one  eddy  cell. 

Let  a flat  plate  airfoil  of  chord  length  2b  be  placed  along  the 
midchord  point  at  the  origin  of  the  x ,y-coordinate  system.  If  6 is 
tion  of  the  £ , n-coordinates  witn  respect  to  the  x ,y-coordinates  and 
tion  of  the  origin  of  the  £ , g-coordinate  system  at  zero  time,  tnen  t 
tional  incoming  flow  to  which  the  airfoil  is  subject  to  is  represent* 
function 


(x-x.-ut)cose  ♦ (y-y--Vt)sine 


Tne  problem  under  consideration  is  to  evaluate  the  unsteady  response  of  tne  flat  plate 
airfoil  subject  to  sucn  an  incoming  flow.  Further  assumptions  made  in  solving  the 
problem  are:  (i)  the  flow  about  the  foil  remains  attached  (ii)  the  introduction  of  the 
foil  into  the  incoming  flow  does  not  change  the  total  circulation  of  tne  whole  system 
and  (iii)  the  foil  remains  stationary. 


STUDY  OF  STREAMFUNCTION 


In  order  to  obtain  an  idea  of  the  streamline  pattern  represented  by  the  assumed 
streamf unction  one  can  assign,  without  losing  any  essential  features,  zero  values  for 
the  orientation  parameters  Xg , yq  and  and  choose  to  plot  the  instantaneous  stream- 
line pattern  at  zero  time.  On  normalizing  the  resulting  streamfunction  with  respect  to 
the  semichord  length  b,  and  the  mean  velocity  components  with  respect  to  the  eddy-inten 
sity  parameter  and  the  semichord  one  gets 


Using  a computer  program  the  stencils  for  the  streamline  patterns  of  ip ' have  been  plotted 
in  Figs.  2 and  3 for  two  sets  of  parametric  values.  Fig.  2 represents  the  streamline 
pattern  for  nontranslating  eddy-array  flow.  For  this  case  the  streamlines  are  also  lines 
of  constant  vorticity. 


It  can  also  be  shown  that  the  assumed  streamfunction:  (i)  satisfies  the  continuity 
equation  for  the  flow  of  an  incompressible  fluid  (ii)  represents  a flow  field  which  is 
in  general  rotational  (iii)  involves  rectangular  eddy  cells  each  of  total  circulation 
strength  equal  to  4g(k/n  + h/k)  (iv)  represents  streamlines  with  2k  and  2h  as  the  wave- 
lengths (the  stencil  sides)  of  the  distance  in  the  x-  and  y-directions  respectively  and 
(v)  the  upwash  created  oy  the  inlet  flow  streamfunction  is  periodic  only  when  the  direc 
tion  of  the  mean  inlet  velocity  satisfies  the  constraint 


cosO 


and  when  6 satisfies  the  above  constraint  the  nondimensional  period  of  the  upwash  is 
given  by 


cose 


PROPOSED  METHOD 


The  method  of  solution  adopted  to  solve  the  present  problem  is  the  singularity  dis- 
tribution principle  in  combination  with  a time-marching  tecnnique.  The  problem  is  solved 
in  two  stages,  namely,  (1)  auxiliary  solution  and  (2)  time-marching  solution.  By  'aux- 
iliary' solution  is  meant  the  solution  of  the  problem  which  completely  neglects  the  pres- 
ence of  the  wake  vortex  sheet  (which  is  always  present  in  such  an  unsteady  flow  situa- 
tion) and  treats  time  t as  a parameter;  hence  in  solving  for  the  auxiliary  solution  the 
unsteady  terms  Ut  and  Vt  in  the  incoming  flow  streamfunction  can  be  lumped  with  the  ori- 
entation parameters  xg  and  yg  respectively.  The  time-marching  part  of  the  solution 
increments  time  by  equal  steps  starting  from  zero  time,  makes  use  of  the  auxiliary  solu- 
tion at  every  time  step,  keeps  track  of  the  shedding  and  growth  of  the  wake  vortex 
sheet,  evaluates  the  unsteady  aerodynamic  responses  (lift  and  momentlat  every  time  step, 
and  continues  along  the  time-axis  up  to  a specified  maximum  time  limit. 

The  question  arises  as  to  whether  one  is  justified  in  solving  the  problem  of  an 
airfoil  located  in  a rotational  flow  field  using  the  method  of  singularities.  Since  it 
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is  assumed  in  the  present  theory  that  the  fluid  is  inviscid  it  is  not  possible  for  the 
airfoil  to  exert  any  nonconservative  force  on  the  fluid  particles  and  thus  change  the 
rotationality  of  the  incoming  flow.  Also,  since  it  is  assumed  that  the  flow  about  tne 
foil  remains  always  attached,  the  only  way  other  than  pure  convective  effects,  by  which 
the  rotationality  of  the  fluid  particles  in  the  flow  field  can  change  is  due  to  the 
shedding  of  wake  vortex  sheet  which  arises  solely  due  to  the  unsteadiness  of  the  total 
bound  circulation  around  the  airfoil.  Since  the  time-marching  technique  employed  takes 
complete  care  of  the  unsteady  wake  vortex  sheet  starting  from  zero  time  the  use  of  con- 
ventional method  of  singularities  (satisfying  the  nonpenetration  and  Kutta  conditions) 
is  justified  in  the  present  work. 

AUXILIARY  SOLUTION 

According  to  the  earlier  stated  definition  of  'auxiliary  problem',  assuming  no  wake 
vortex  sheet  to  be  present  and  considering  time  t only  as  a parameter,  the  upwash  along 
tne  airfoil  is  to  be  calculated  from  the  streamf unction 

Ily(x,y)  » Uy  - Vx  ♦ g cos£  [(x-xQ)cose  ♦ (y-ygisine]  cos£-  [-(x-x0)sin6  ♦ (y-yg)cose] 


where  suffix  A denotes  'auxiliary';  Xq  and  yg  incorporate  in  them  the  temporal  terms  Ut 
and  Vt  also.  For  the  sake  of  simplicity  it  is  chosen  not  to  indicate  by  symbols  in  the 
foregoing  and  following  expressions  the  fact  that  time  t is  a parameter. 

The  upwash  (normal  velocity  component)  along  the  foil  becomes 


vA(x,0) 


8x  y*0 


* V ♦ *£cose  sinj-  [(x-x0)cos8  - y0sine]  cos£-  [-(x-x0)sine  - ygCose] 

- gjUine  cos£  [(x-Xq) cosB  - y0sine]  sin£  [-(x-x0)sin8  - yQcose] 

Decomposing  the  upwash  into  one  due  to  steady  translation  (denoted  by  suffix  S)  and  one 
due  to  the  presence  of  eddies  (denoted  by  suffix  E)  as  vA  = Vg  + Vg  one  gets  after  con- 
siderable simplification 


vs(x,0)  = V 

vL(x,0)  « ^ ^Z^cosZjSin(Z2  £■)  - ZjSinZjCos^  + ZjCOsZ^sinCZj  jp)  - ZjSinZ^cosCZj  g-)J 


The  new  constants  involved  in  the  above  equation  have  the  following  expressions 

,b b 


Zj  « cos0  » — sin8) 


Z2  * *(j-  cose  - jj-  sine) 


Zj  * ir(^  cos 8 ♦ sine) 


■ ■ [?]> 


b b 

z4  * *(jj-  cosQ  - jj-  sin0) 


z6  * If0]2!  ■ [?h 


I 


tion 


Using  the  Sohngen  inversion  formula  (Ref.  1)  to  solve  the  singular  integral  equa- 


! 


1 


cos* 


where  the  variable  * has  been  defined  by  the  expression  x : b cos*.  It  may  be  noted  that 
even  tnough  the  expression  for  the  auxiliary  bound  vortex  distribution  involves  infinite 
series  its  zeroth  and  higher  moments  about  the  midchord  have  expressions  wnich  involve 
only  finite  number  of  terms.  Only  these  moments  are  needed  in  the  time-marching  solu- 
tion. 


TIME -MARCHING  SOLUTION 


In  the  present  analysis  the  unsteadiness  of  the  resultant  flow  field  has  three 
gins:  (i)  unsteady  nature  of  the  incoming  flow;  this  is  taken  care  of  by  the  stream 
tion  assumed  (ii)  unsteadiness  of  the  bound  vortex  distribution  and  (iii)  unsteady 
ture  of  the  shedding  and  growth  of  the  wake  vortex  sheet.  In  order  to  take  care  of 
unsteadiness  of  the  bound  vortex  distribution  and  the  wake  vorte;;  sheet  a time-marc- 
scheme  can  be  devised  as  follows . 


decomposing  (Ref.  2,  3)  the  bound  vortex  distribution  and  the  total  circulation  in 
to  two  components,  namely,  one  due  to  'auxiliary'  solution  (solution  of  flow  entirely 
neglecting  the  presence  of  wake)  denoted  by  suffix  A and  another  due  to  'unsteady'  solu 
tion  (solution  of  flow  entirely  due  to  the  presence  of  wake)  denoted  by  suffix  U 


Making  use  of  the  condition  that  the  presence  of  the  foil  (and  hence  the  presence  of  the 
wake)  do  not  change  the  total  circulation  of  the  whole  system  one  gets 


where  W(t)  denotes  the  entire  instantaneous  line  element  representing  the  wake  vortex 
sheet;  S is  the  arc  length  along  the  line  element  W(t)  measured  from  the  trailing  edge 
of  the  foil;  and  yw(S)  is  the  wake  vortex  distribution  (Fig.  4). 

Rearranging  Eq.  (1)  and  using  Eq.  (2) 


where  K(§,x)  denotes  the  bound  vortex  distribution  along  tne  foil  due  to  the  presence 
of  a unit  vortex  at  vector  distance  S from  the  trailing  edge;  S is  the  vector  from 
trailing  edge  associated  with  a point  on  the  vortex  sheet  whose  arc  length  is  S (Fig. 
Changing  the  order  of  integration  and  grouping  terms 


means  the  total  circulation 


after  defining  a new  kernel  to  be  G- 

around  the  foil  due  to  a unit  vortex  located  at  vector  distance  S 
sions  for  K(S,x)  and  Gs  are  derived  in  Appendix  A. 


The  above  integral  equation  for  wake  vortex  strength  will  be  solved  in  a time- 
marching fashion  making  use  of  (i)  the  concept  (Ref.  3)  that  the  wake  which  is  a con 
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tinuous  distribution  of  point  vortices  may  be  approximated  by  an  equivalent  discrete 
distribution  of  point  vortices;  the  discretization  scheme  adopted  in  this  work  is  to 
represent  the  elemental  wake  vortex  sheet  shed  during  any  one  time  step  by  a single 
point  vortex  of  strength  equal  to  the  total  circulation  around  the  element  and  posi- 
tioned at  a meanpoint  of  the  element;  (ii)  the  condition  that  the  wake  vortices  are 
free  in  the  sense  that  they  move  with  the  fluid  element  to  which  they  are  'attached'. 

The  details  of  the  solution  scheme  follows. 

At  time  t+At  denoting  the  far  end  (Fig.  4)  of  the  wake  vortex  sheet  (i.e.,  location 
at  time  t+At  of  the  fluid  particle  which  marked  the  birth  of  the  wake  vortex  sheet  at 
time  t = 0)  by  Sglt+At)  one  gets  from  the  wake  vortex  integral  equation 

b 

% Mi  * Gsi  vs> ds 


Measuring  time  in  steps  of  equal  length  At  and  dividing  the  wake  vortex  sheet  into  seg- 
ments formed  during  each  time  step  (Fig.  4) 


rA(n*l  At) 


n Sj (n+lAt ) 

21  / [l  + G ]y  (S)  dS 

1 = 1 Sl-l(n-7lAt) 


/ [1  * G ]y  (S)  dS 
S CrTTlAt)  s ,v 


where  S^tn+l  At)  denotes  the  location  at  time  n+1  At  of  the  fluid  particle  which  marked 
the  birth  of  the  wake  vortex  sheet  at  time  lAt. 


Discretizing  the  wake  segments,  that  is,  replacing  each  segment  of  continuous  vor- 
tices by  a single  vortex 

-rA(n^At)=  |j[l+G1s]g1  + [l+G-1]^ 


where  g^  denotes  the  constant  circulation  strength  (equal  to  the  total  vorticity  of  the 
fluid  particles  constituting  the  elemental  wake  segment  formed  during  the  time  interval 
l-lAt<t<lAt)  of  the  1th  youngest  discrete  wake  vortex  which  is  located  at  a mean  point 
of  the  wake  segment  it  represents.  Thus 

-rA(STI  At)  - [l  ♦ G*],, 

*n+1  TT^f1 


where  represents  Gs  when  the  arc  length  involved  in  the  evaluation  of  Gg  is  the  one 
with  respect  to  the  location  of  the  1th  youngest  wake  vortex. 

The  location  vector  Sj  (n+1  At)  for  the  1th  youngest  wake  vortex  at  time  iv+T  At  may 
be  found  by  the  following  predictor-corrector  scheme  adoDted  from  Ref.  3 


Pj  = §j(nAt)  ♦ At  Vj (§j (nAt) ) 

„ _ „ fv.  ($.  (nAt))  ♦ V.  (P  ) 

Sj(n*lAt)  = Sj(nAt)  ♦ At  * 1 1 1 


where  V}  is  the  convective  velocity  for  the  1th  wake  vortex.  Obviously,  such  a predictor- 
corrector  scheme  is  necessary  due  to  the  fact  that  the  val_j  of  the  convective  velocity 
keeps  changing  with  the  convection.  The  convective  velocity  v^  for  the  1th  wake  vortex 
is  composed  of  the  following  components 


f' 


1 


V 


freestream 


^auxiliary  bound 
vortices 


unsteady  bound 
vortices 


where  the  summation  index  k stands  for  all  wake  vortices  except  the  1th.  The  term  ' free- 
stream'  includes  the  steady  translation  and  eddy-array  components.  Assuming  that  the 
elemental  portion  of  the  wake  vortex  sheet  being  shed  does  not  contribute  to  the  con- 
vection of  the  wake  all  the  information  needed  to  calculate  Vj^  is  known.  In  the  limit 
of  small  time  steps  this  assumption  is  valid. 

Due  to  the  limitation  of  the  computer  memory  it  is  not  possible  to  represent  the 
wake  vortex  by  more  than  a reasonably  large  number  of  discrete  free  vortices.  This 
calls  for  a scheme  by  which  one  will  keep  track  of  only  a finite  number  of  wake  vortices 
and  allow  for  a rational  "overflow"  of  the  remaining  wake  vortices.  Such  a scheme  fol- 
lows. It  may  be  seen  from  the  results  of  Appendix  A that  for  large  values  of  the  loca- 
tion vector  3 the  kernel  Gg  approaches  zero.  Hence  if  it  is  agreed  to  handle  a maximum 
of,  say,  only  N number  of  wake  vortices  and  it  is  expected  that,  in  general,  the  older 
wake  vortices  would  have  been  swept  away  farther  from  the  foil  the  time-marching  equa- 
tion for  the  strength  of  (n+l)th  nascent  wake  vortex  becomes 


UNSTEADY  RESPONSE  EXPRESSIONS 


Starting  from  Euler's  equation  of  motion  it  can  be  shown  that  for  the  flow  of  an 
incompressible,  inviscid  fluid  about  a flat  plate  airfoil  the  loading,  lift  (positive 
upward) , and  moment  about  the  leading  edge  (positive  clockwise)  are  given  by 


/ (j  <-bx)y(x,t)  dx 


where  u(x,t)  is  the  resultant  tangential  velocity  along  the  foil  position  and  y(x,t)  is 
the  bound  vortex  distribution. 

Decomposing  the  tangential  velocity  and  bound  vortex  distribution  into  steady  trans 
lation  (suffix  S) , eddy-presence  (suffix  E)  and  unsteady  (suffix  U)  parts  and  grouping 
the  response  expressions  into  auxiliary  (suffix  A) , unsteady  (suffix  U)  and  mixed  (suf- 
fix M)  contributions  one  gets  after  omitting  to  indicate  the  dependency  on  independent 
variables 


U / (x*b)yF  dx  - / (x+b)UpY_  dx 


nent  coefficients  to  be 


one  gets  the  different  contributions  to  be 


where 
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C1A  * C1A1  * C1A2  * C1A3  * C1A4  * C1A5  * C1A6 

1 111  * C1U1  * C1U2  * C1U3 

C1M  * C1M1  * C1M2  * C1M3  * C1M4 

^rnA  ~ ^mAl  * ^“mA2  * ^mA3  * ^mA4  "mA5  ^mA6 

l'mU  * CmUl  * SnU2  * SnU3 

CmM  * Ci*ll  * CmM2  * CmM3  * CmM4 


Nondimensionalizing  the  time  t,  distance  x,  bound  vortex  strength  y,  and  the  tan- 
gential velocity  u according  to  the  following  definitions 


t' 


' b • 


ub 

8 


one  gets  after  arranging  each  of  the  expressions  involved  to  be  the  product  of  nondimen- 
sionalized  terms 


Sai  - f [fw]2  / YS  - 


C1U2  * g [bwj  at'  [ Yu  dx' 


"1A2  LbW. 
Ub 


hwl  / UEYE  dx 


-1 


Ub  ^ l2  3 , J1*.,.  dx’ 


1 


“1U3  g Lbw  J at 


-l 


1A3 


'1A4 


IAS 


g I5wj  at’  ye  dx’ 
, i 


Ub 


Ub 


, / x’y’  dx’ 
-1  c 


C1M1  “ [bi*]  | dx’ 

2 1 

C1M2  * [bw]  ( UUYE  dx’ 


/ T p djl’ 
-1  fc 


C1A6  ’ [bw]  .{  UE  YS  dX' 

2 1 

CIUI  ■ [fo]  _{  “u  Yu  dx’ 

CmAl  * ' }f  Bw]  [.{  *'YS  dX'  * J YS  dx’] 


-1 


Ub  g 1 r , . , 
"1M3  g [bwj  YU  dx 


C1M4  * §w]  uEYU  dX' 


2 1 


CmA2  * - Mbw]  [ [ x'uEyE  dx'  * J uEyE  dx'l 

CmA3*-|r[bw]  [ It’  _{  yE  dx'] 

CmA4  =?jT  [bwj  [j  at"'  _{  x'  ye  dx'  * at'  | X'yE  dx'] 
CmA5  * - irlSw]  [_{  X’yE  dx'  * _{  yE  dx'] 

CmAft  ’ ' 2 [bw]  [ X'UEyS  dx ' UEYS  dx '] 

CmUl  * - KGw]  [ _{  X'%YU  dx'  * _{  UUYU  dx'] 

c'mu2 1 - ? riGw]  [k-  \ YU  dx'] 

cmU3  ’ 2 “gfbW ] [lat'j  x'  Yu  dx'  * IT'  X'YU  dx'l 


c*mi  * - l[bw]2  [ A'Vs  dx'  * .{  “uYs  dx’J 

S*I2  * - l[bw]2  [_{We  dX’  * .{Ve  dX’^ 

CmM3  - ’ ? dX’  *_{  YU  H 

CrrW4  * ' Mivf]2  dX’  + J UEYi)  dX’l 


Using  the  results  of  the  auxiliary  solution  for  the  bound  vortex  distributions  and 
the  following  expressions  for  the  tangential  velocity  components 


Ug  = h [Z^cosZ^sintZjX')  - Z^sinZ^cosfZjX')  - ZjCosZjSinfZjX’)  ♦ ZjSinZjCos (ZjX')] 


one  gets  after  carrying  out  the  integrations  involved 


/ y'  dx’  = 2v 


/ x'y'  dx'  = -ir  — 


/ y£  dx'  * nfZ^jCZjJcosZj  ♦ ZjJj  [Z1)cosZ6  - Z2JQ (Z2)sinZg  - ZjJ0(Zj)sinZ6] 

1 

/ x’Yg  dx'  = x[J j (Z2)sinZj  ♦ Jj(Zj)sinZ^] 

/ X'2  y’  dx'  = \ [t(Z2-  ^)J1(Z2)  ♦ 2J0(Z2))cosZ5  ♦ { CZr  j^JjCZj)  ♦ 2-yZ,,) }cosZ6 
- ZjJpCZjisinZj  - ZjJ^ZjlsinZ^l 

/ Yj  dx'  * j [f  Z2Jfl(Z2)sinZ5  + Z1JQ(Z1)sinZ6)  { Z^osZ^inZj  - ZjCOsZjSinZj) 

♦ 3 { Z2JjfZ2)cosZ5  ♦ ZjJl(Z1)cosZ6)  { JgCZjJZjSinZj  - J^Z^sinZ^ 

♦ { Z2Jj(Z2)cosZ5  ♦ ZjJj(Zj)cosZg)  { Z4sinZ6cosZ1  - ZjSinZjCosZj) 

♦ 2 ( Z2J0CZ2)sinZs  ♦ ZJJ0(ZJ)sinZ6)  { J0(Zj)Z4sinZ6  - JpfZjJZjSinZj}] 

1 

/ x'UgYg  dx'  = g-  1 2 { Z2JQ(Z2)sinZ5  ♦ ZjJ0(Zj)sinZ6)  { Z3cosZjSinZ2  - Z^osZgSinZj } 

* (Z2J0(Z2)sinZs  ♦ Z1J0(Z1)sinZ6l  {Jgtf^ZjSinZj  - .yZjJZ^inZ^ 

♦ ( Z2JQ(Z2)slnZj  ♦ ZjJ()(Z1)sinZ6)  { ZjSinZjCOSZj  - Z^inZ^osZj)  ] 

/ UgVj  dx’  * y f-Z4cosZ6sinZj  - 2J0(Zj)Z4slnZ6  ♦ ZjCOsZjSinZj  ♦ 2J0 (ZjJZjSinZjJ 

/ x'UgYg  dx'  = J I2Z4cosZgSinZj  ♦ Z4*inZ6  {Jq(Zj)  ♦ cosZj)  - ZZ^cosZjSinZj 
- ZjSinZj  ( JQ(Z2)  ♦ cosZj}  ] 

Other  integrals  involved  in  the  expressions  for  the  various  components  of  the  re- 
sponse coefficients  cannot  be  evaluated  in  closed  form;  they  have  to  be  obtained  by 
'■••tun.  This  calls  for  a special  scheme  to  handle  in  closed  form  the  integration 


.-»•  - ■' 


I 


30-10 


over  a small  region  around  the  leading  edge  which  is  a singular  point.  Adopting  such  a 
scheme  based  on  the  fact  that  close  to  the  leading  edge  (x'  = -1)  one  can  approximate 
YU  by  C/  /x'+l  where  C is  a constant  and  making  use  of  the  following  expression  for  the 
unsteady  component  of  the  tangential  velocity  along  the  foil  position 


% 


.-.-2 


n = Wake 

vortices 


g ' 

7T 


- x')2  * Xu 


one  gets 


1 1 
J Yy  dx'  =<  2e  YyC-l+e)  ♦ / Y u dx' 


-1+c 


7 1 

J x'Yy  dx'  = - j e (3-e)  Yy(-l+e)  + / x'Yy 


dx' 


/x'^  dx'  = ^ eUS-10e+3e2)  Y^C-l+e)  + J x'2Y[)  dx- 


1 

/ dx-  - 

- b X 

n 

f x'u^  dx- 

1 

2t 

n 

{ “uYs  dx'  * 

vb  y _ 

g n x’: 
n 

1 

/ x'UyYg  dx' 

1 Vb 

2 g Z_  ' 

6 n j 

J UyYg  dx'  - 

i 

’ 2x 

X- 

n x 

i 

1 

J x'u^  dx- 

1 

’ ‘ Ti 

[?; 

2e  ..j1 


Yu 


dx1 


(x'+l)2  ♦ y'2  -l*e(x'-x')2  ♦ y'2 


- (2/3)e  (3-e)  YU(~1*t:) 

(x’+l)2  + y'2 
v n J Jn 


♦ s — 

-l*e  (x 


dx' 


p 

yn  J 


J Yp  dx' 


-1 


g'y'  1 
“n'n  1 

,2 


/ x'yi  dx’ 


-1 


1 ..I 

/ UgYy  dx'  =>  EY^(-l+t)  [{ -Z4sin(Z1+Z6)  + ZjSinCZj+Zj)}]*  / u^Yy  dx' 

J x ' UgYy  dx'  =>  - j t(3-c)  Y,)(-l+e)[{ -Z4sin(ZJ+Z6)  ♦ zjsin(Z2*V^*  \ x’u^Yy  dx' 


‘ 


where  e is  a small  number  very  much  smaller  than  unity. 

LIMITING  CASES 

In  the  present  theory  if  the  eddy- intensity  parameter  is  made  zero  it  is  possible 
to  get  the  classical  result  for  the  case  of  a flat  plate  airfoil  in  a uniform  parallel 
flow.  By  assuming  the  length  of  one  side  of  the  rectangular  eddy  to  be  very  large  the 
present  work  can  simulate  stratified  flows  over  thin  foils.  It  can  be  seen  from  the 
assumed  inlet  flow  streamfunction  that  the  freestream  velocity  components  along  the  foil 
position  for  the  case  of  0=0,  xq=0,  yg=0  and  V=0  becomes 

v(x,0,t)  = g£  sin£(x-Ut) 
u(x,0,t)  = U 

Thus  for  this  special  set  of  parametric  values  the  upwash  is  that  of  a sinusoidal  gust 
and  if  the  quantity  gir/k  is  relatively  small  in  comparision  with  U then  the  formulation 
is  exactly  the  same  as  that  of  the  classical  transverse  sinusoidal  gust  problem.  How- 
ever it  should  be  noted  that  the  present  theory  does  not  (i)  limit  itself  to  a perturba- 
tion type  of  analysis  (ii)  assume  the  spatial  location  and  temporal  variation  of  the 
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wake  vortex  sheet  (iii)  assume  the  temporal  variation  of  the  bound  vortex  distribution; 
these  are  some  of  the  modeling  assumptions  in  the  existing  theories  for  sinusoidal  gust 
problems.  Another  important  difference  to  be  noted  is  that  in  the  sinusoidal  gust  the- 
ories nothing  is  said  explicitly  about  the  velocity  distribution  of  the  freestream  flow 
field  except  the  specification  of  the  upwash  whereas  in  the  equivalent  limiting  case  of 
the  present  work  the  flow  field  continues  to  be  that  of  a translating  rectangular  grid 
of  nonpotential  eddies.  Hence  it  is  not  expected  that  the  results  of  this  limiting  case 
will  be  identical  to  those  of  the  sinusoidal  gust  theories  even  after  allowing  for  the 
modeling  simplifications  noted. 

RESULTS  AND  RECOMMENDATION 

A computer  program  has  been  developed  to  obtain  numerical  results  for  the  response 
coefficients  making  use  of  the  expressions  derived  in  this  theory.  The  results  of  this 
program  is  compared  with  the  linear  theory  of  Sears  (Ref.  6)  in  Table  1.  k is  the  re- 
duced frequency  and  vq  is  the  gust  amplitude.  In  Table  2 the  effect  of  large  amplitudes 
of  sinusoidal  gust  on  the  response  is  given.  It  has  been  checked  that  the  program  gives 
satisfactory  values  for  uniform  parallel  flow  response  when  the  eddy  strength  is  made 
very  small.  Figure  6 shows  the  shape  of  trailing  vortex  sheets  at  the  instant  when  the 
foil  has  traveled  through  four  sinusoidal  gust  wavelengths.  It  is  intended  to  develop  a 
similar  analysis  for  the  case  of  a cascade. 


APPENDIX  A 

Consider  a flat  plate  airfoil  of  chord  length  2b  to  be  located  at  the  origin  (Fig. 
5)  of  Z-plane  along  the  real  axis  in  the  flow  field  of  a vortex  of  strength  P located 
at  a general  point  Zq.  The  expressions  for  the  bound  vortex  distribution  and  the  total 
circulation,  subject  to  nonpenetration  and  Kutta  conditions,  may  be  derived  as  follows. 
The  conformal  mapping  function 


maps  the  flat  plate  into  a circle  of  radius  b as  shown  in  Fig.  5.  Let  £q  be  ihe  m 
of  Zq . Consider  a vortex  of  strength  -P  to  be  placed  at  the  'image'  point  b2/7o. 
complex  potential  for  the  system  in  the  transformed  plane  is 

v2 


mapping 
The 


v(C)  ■ i 2^  log (C-Cn)  - i 57  logCt- 


2tt 


) 


Making  use  of  the  following  expression  for  the  complex  conjugate  velocity 

dw 

u - 1V  * 37 

and  the  fact  that  on  the  circle  under  consideration  £ = b(cos0  + i sine)  the  expressions 
for  the  Cartesian  components  of  velocity  can  be  shown  to  be 


Id  *b 


r 

2it  b 


tei 


2 2 

n„  - b ) sine 


2 2 2 

K0  ♦ 13  * b - 25 gb  cose  - 2rigb  sine) 


\l  *b 


r_  i 

2n  b ,„2 


k; 


2 i2-. 

n„  - b ) cose 


2 2 

C5‘  ♦ n0  ♦ b - 25Qb  cose  - 2nQb  sine) 


by 


Measuring  9 anticlockwise  from  the  £-axis  the  polar  components  of  velocity  are  given 


Vf  » u cose  + v sine 
V0  = -u  sine  + v cose 


and  hence 


9 * 2nb 


r2  2 l2 

so  * "o  * b 


(5 


q - n0  * b - 25gb  cose  - 2ngb  sine) 


This  proves  that  the  circle  |c|=b  is  a streamline  and  hence  the  satisfaction  of  the  non- 
penetration condition. 

At  the  point  9-0  corresponding  to  the  'trailing  edge'  Z=b  the  tangential  velocity  is 


le*o 


r 

2^77 


,2  2 >.2 

s0  * n0  ‘ b 


50  * nQ  * b - 25 Qb 


■a—-.. * 
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Adding  a circulation  to  make  the  tangential  velocity  vanish  at  the  trailing  edge  (Kutta 
condition)  and  resubstituting  b cos8=£;  b sin9=o 


V = X 
6 2 rib 


.2  2 1.2 
Sn  + % - b 


s0  * - b‘ 


to  * V h2  - ■ 2V  «0  * \ * 1)2  - 25Ob 


From  page  46  of  Ref.  5 one  gets  the  relation  between  the  bound  vortex  distribution  and 
the  tangential  velocity  component,  after  noting  the  difference  in  the  transformation 
functions  used,  to  be 


Y(€;n) 


sin9 


2bV 


6 


Hence,  using  the  relation  n = / b2-C2  for  the  upper  half  of  the  circle,  one  gets  for  b=l 


YCO  = W 


r'  «o  * bo  * 1 


7TT 


i 


i 


C2  - ho  ♦ 1 ' 25o  ’ S**  n2  ♦ i - 2So5  - V1  - 


The  kernel  function  K(S,x) needed  in  the  time-marching  solution  takes  the  form 

2 

no 


* 1 * % - 1 
K(S,x)  = I 


1 1 

L?0  * no  * 1 ■ 25o  S2  ♦ ^ i - 2V  - ^ ^ 


where  S corresponds  to  xg  iyg  and  by  inverting  the  mapping  function  one  gets 


S * x 

S0  * x0  * cosf 


n0  = y0  ♦ /R  Sin| 


R -M  ' yl  - * 4x‘y, 


0 * Tan 


-1 


)2 

Vo 


2 2 
0y0 


2 2 , 

xo  - yo  - 1 


where  Tan'2'  stands  for  principal  value. 

For  the  particular  case  of  no=0  the  result  for  the  bound  vortex  distribution  re- 
duces to 


Y(C) 


«0  * 1 
o 2V  5°  ' 1 


- £’  I1  - s 2Co 

' n jTTl?  , + ,2  _ 

_ i"  /i  - c i /xo  * 1 
■ it  v i * e r-^r 


- s 1 

5 x0  - -0 

This  agrees  with  the  result  reported  in  Ref.  2. 

The  total  circulation  around  the  flat  plate  foil  for  b=l  is  given  by 


I 

[! 

i 


r«0,n0)  * / y(C)  d? 


♦ bo  ' 1 


l5o  ♦ *0  + 1 ' 25o 


5o  * no  ' 1 


Hence  the  kernel  function  Gg  needed  in  the  time-marching  solution  becomes 


jLxLj_ 

?o  * ”0  * 1 ' 25o 


- 1 - 2 Tan'1 

TT 


2rio 

Zl  2 T 

eo  * no  • 1 


For  the  particular  case  of  no=0  one  gets  from  above  result  for  total  circulation 

.[s0  + 1 

ra  ) * r — 1 

which  also  entirely  agrees  with  the  result  of  Ref.  2. 


w 
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Table  I 

Sears  Function  S and  Fundamental  of  Present  Result  P for  vQ/U 


= .01 


k 

|S(k) | 

ife(k) 

1 P (k)  | 

/P(k) 

0.02 

0.96583 

-4.439 

0.96498 

-1.906 

0.05 

0.91431 

-8.060 

0.91867 

-4.289 

0.10 

0.83735 

-11.258 

0.85440 

-7.368 

0.50 

0.52647 

-4.797 

0.56922 

-19.385 

1.00 

0.38957 

18.861 

0.29153 

-32.947 

5.00 

0.17820 

-117.097 

0.69366 

2.331 

Table  II 

Effect  of  Gust  Amplitude  on  Response 


T^7UT 


Ao  + 

+ 

+ 


A^cos  ut  + 
AjCOs  5iut 
B^sin  3wt 


A2cos  2wt 
+ AgCos  6ut 
+ B4sin  4wt 


+ A^cos  3 dit  + A^cos  4cut 
+ B^sin  wt  + BjSin  2u)t 
+ BjSin  5wt 


D 

S 

o 

> 

onn 

o.oS 

571(5 — 

0.50 

1.6(5 

A„ 

-0.00095 

0.0031$ 

5755S75- 

0.00731 

0.00383 

A? 

0.68843 

0.68861 

0.68932 

0.70759 

0.73617 

A* 

-0.00090 

-0.00400 

-0.00783 

-0.03637 

-0.06334 

A? 

-0.00013 

-0.00041 

-0.00116 

-0.02216 

-0.07056 

A l 

-0.00012 

-0.00017 

-0.00019 

0.00083 

0.00613 

4 

-0.00011 

-0.00019 

-0.00030 

-0.00734 

-0.05103 

4 

-0.00011 

-0.00020 

-0.00036 

-0.01065 

-0.07169 

B , 

-5.32406 

-5.32486 

-5.32668 

-5.35415 

-5.38706 

-0.00063 

-0.00254 

-0.00488 

-0.02338 

-0.04263 

4 

-0.00008 

0.00016 

0.00107 

0.03087 

0.12100 

4 

-0.00005 

-0.00013 

-0.00029 

-0.00441 

-0.01860 

4 

-0.00001 

-0.00009 

-0.00026 

-0.00534 

-0.02160 

m 


TECHNICAL  REPORT  DISTRIBUTION  LIST 

Stevens  Institute  of  Technology  Contract  NO 0 01 4- 6 7 -A- 020 2- 0016  - NR  094-363 


Turbomachine  Blade  Flutter 


Recipient 


No.  of  Copies 


Defense  Documentation  Center,  Bldg.  5 
Cameron  Station 
Alexandria,  Virginia  22314 


Chief  of  Naval  Research 
Department  of  the  Navy 
Arlington,  Virginia  22217 

Attn:  Mr.  J.  R.  Patton,  Jr.,  Code  473 


Commander 

Naval  Air  Systems  Command 
Department  of  the  Navy 
Washington,  D.C.  20361 

Attn:  Mr.  R.  R.  Brown,  AIR  330 

Dr.  H.  Rosenwasser,  Code  31 0C 
Technical  Library  Division,  Code  604 


Commander 

Naval  Ship  Systems  Command 
Washington,  D.C.  20360 

Attn:  Mr.  R.  R.  Peterson,  Code  03413 

Mr.  C.  Miller,  Code  6146 
Technical  Library 


General  Electric  Company 
Aircraft  Gas  Turbine  Division 
Cincinnati,  Ohio  45215 

Attn:  Manager  of  Engineering 


Commanding  Officer 
U.S.  Army  Research  Office 
Box  CM,  Duke  Station 
Durham,  North  Carolina  27706 
Attn:  ORDOR-PC 


Westinghouse  Electric  Corporation 

Steam  Turbine  Division 

Lester  Branch,  Philadelphia,  PA. 


(1) 


I 

1 1 


United  Aircraft  Corporation 

Pratt  and  Whitney  Aircraft  Division 

P.0.  Box  2691 

West  Palm  Beach,  Florida  33402 


Varo,  Incorporated 
402  East  Gutierrez  Street 
Santa  Barbara,  California  93101 
Attn:  K.J.  Widiner 


Director 

Ordnance  Research  Laboratory 
Pennsylvania  State  University 

University  Park,  Pennsylvania  16802  (1) 


Commander 

Naval  Ordnance  Systems  Command 
Department  of  the  Navy 
Washington,  D.C.  20360 

Attn:  Mr.  B.  Drimmer,  Code  033 


Naval  Ship  Research  and  Development  Center 
Annapolis  Division 
Annapolis,  Maryland  21402 
Attn:  Library,  Code  A214 


U.S.  Naval  Oceanographic  Office 
Suitland,  Maryland  20390 
Attn:  Library,  Code  1640 


Director 

National  Aeronautics  and  Space  Administration 

Headquarters 

Washington,  D.C.  20546 

Attn:  Division  Research  Information 


(1) 


(1) 


(1) 


(1) 


Director 

National  Bureau  of  Standards 
Gaithersburg,  Maryland  20760 


Office  of  the  Assistant  Secretary  of  Defense 
Room  3E1065  - The  Pentagon 
Washington,  D.C.  20301 
Attn:  Technical  Library 


Aerojet-General  Corporation 
P.0.  Box  296 

Azusa,  California  91702 
Attn:  Librarian 


AiResearch  Manufacturing  Company 
9851  Sepulveda  Boulevard 
Los  Angeles,  California  90045 
Attn:  Chief  Engineer 


General  Motors  Corporation 
Allison  Division 
Indianapolis,  Indiana  46206 
Attn:  Director  of  Engineering 


General  Electric  Company 
Engineering  Department 
Turbine  Division 
Schenectady,  New  York  12305 


General  Electric  Company 
Aircraft  Turbine  Department 
West  Lynn,  Massachusetts  01905 


AVCO  Corporation 
Lycoming  Spencer  Division 
652  Oliver  Street 
Williamsport,  Pennsylvania 


17701 


National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  Ohio  44135 


United  Technologies  Corporation 
Pratt  and  Whitney  Aircraft  Division 
East  Hartford,  Connecticut  06118 
Attn:  Chief  Engineer 


Solar  Aircraft  Company 
San  Diego,  California  92101 
Attn:  Chief  Engineer 


SMVM 


Recipient 


No.  of  Copies 


United  Technologies  Corporation 
Research  Department 
East  Hartford,  Connecticut  06118 
Attn:  Director  of  Research 


Commander 

U.S.  Naval  Ordnance  Laboratory 
White  Oak 

Silver  Spring,  Maryland  20910 
Attn:  Library 


Commander 

Naval  Weapons  Center 
China  Lake,  California  93555 
Attn:  Technical  Library 


Commander 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base,  Ohio  45433 
Attn:  AFAPL/RJ 


Director 

U.S.  Naval  Research  Laboratory 
4555  Overlook  Avenue,  S.W. 

Washington,  D.C.  20390 

Attn:  Technical  Information  Division 

Library,  Code  2029  (ONRL) 


Director 

Office  of  Naval  Research  Branch  Office 
1030  East  Green  Street 
Pasadena,  California  91101 


Officer  in  Charge 
Naval  Ship  Engineering  Center 
Philadelphia  Division 
Philadelphia,  Pennsylvania  19112 
Attn:  Code  6700 

Technical  Library 


Bureau  of  Naval  Personnel 
Department  of  the  Navy 
Washington,  D.C.  20370 


Recipient 


No.  of  Copies 


Navy  Underwater  Sound  Laboratory 
Fort  Trumbull 

New  London,  Connecticut  06320 

Attn:  Technical  Library  (1) 


U.S.  Naval  Weapons  Laboratory 
Dahlgren,  Virginia  22448 

Attn:  Technical  Library  (1) 


Director,  Project  SQUID 

Jet  Propulsion  Center 

School  of  Mechanical  Engineering 

Purdue  University 

Lafayette,  Indiana  47907 

Attn:  Dr.  S.N.  Murthy  (1) 

Dr.  Bruce  Reese  (1) 


Army  Missile  Command 

Research  and  Development  Directorate 

Redstone  Arsenal,  Alabama  35809 

Attn:  Propulsion  Laboratory  (1) 


Director 

Applied  Physics  Laboratory 
8621  Georgina  Avenue 
Silver  Spring,  Maryland  20910 

Attn:  Library  (1) 


Commander 

U.S.  Air  Force  Systems  Command 
Andrews  Air  Force  Base 
Silver  Hill,  Maryland  20331 


(1) 


Commander 

Air  Force  Aero  Propulsion  Laboratory,  TBP 
Wr ight-Patterson  Air  Force  Base 

Dayton,  Ohio  45433  (1) 


Commander 

Air  Force  Rocket  Propulsion  Laboratory 
Edwards  Air  Force  Base,  California  93523 

Attn:  AFRPL  (2) 


Recipient 


No.  of  Copies 


Naval  Missile  Center 

Point  Mugu,  California  93041 

Attn:  Technical  Library  (1) 

Code  5632.2 


Headquarters 
Naval  Material  Command 
Special  Projects  Office 
Washington,  D.C.  20360 

Attn:  Technical  Library  (1) 


Commander 

Air  Force  Office  of  Scientific  Research 
1400  Wilson  Boulevard 
Arlington,  Virginia  22209 

Attn:  J.F.  Masi  (1) 


Directof  of  Defense  Research  and  Engineering 

Technical  Library 

Room  3C128  - The  Pentagon 

Washington,  D.C.  20301 

Attn:  Propulsion  Technology  (1) 


Chief  of  Research  and  Development 
Headquarters,  Department  of  the  Army 
Washington,  D.C.  20310 
Attn:  Dr.  S.  J.  Magram 

Physical  & Engrg.  Division  (1) 


Dr.  Larry  D.  Mitchell 
Associate  Professor 

Department  of  Mechanical  Engineering 
Virginia  Polytechnic  Institute  and 
State  University 
Randolph  Hall 

Blacksburg,  Virginia  24061  (1) 


Dr.  R.  H.  Kemp 
NASA-Lev/is  Research  Center 
MS  49-3 

21000  Brookpark  Road 
Cleveland,  Ohio  44135 


(1) 


' 11 


MWMIMRai 


Recipient 


Mo.  of  Copies 


Assistant  Chief  for  Technology 
Office  of  Naval  Research,  Code  200 

Arlington,  Virginia  22217  (1) 


AVCO  Lycoming  Division 
550  South  Main  Street 
Stratford,  Connecticut  06497 

Attn:  Division  Library  (1) 


Mr.  F.  E.  Walters 

Office  of  Naval  Research  Branch  Office 
715  Broadway  (5th  Floor) 

Hew  York,  N.Y.  10003  (1) 


i 


security  classification  of  this  page  (Mien  D«t«  Entered) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


REPORT  DOCUMENTATION  PAGE 


3.  RECIPIENT'S  catalog  number 


t.  REPORT  NUMBER 


2.  GOVT  ACCESSION  NO. 


5.  TYPE  OF  REPORT  4 PERIOD  COVERED 


TITLE  (and  Subtitle) 


Research  on  the  Flutter  of  Axial 
Turbomachine  Blading*,  -- 


Reprint  1976 


*■  PERFORWIYtJ" 


ME-RT-760fQ/3  ; - 

-fF  contract  or  graht  nuT’ be r<c«; 


7.  AUTHORfsJ 


N0dOli-76 


Sisto7  P.V.K.  Perumal 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


AREA  4 WORK  UNIT  NUMBERS 


Stevens  Institute  of  Technology 
Castle  Point  Station 
Hoboken,  N.J.  07030 


II.  CONTROLLING  OF  FICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 
715  Broadway  (5th  Floor) 

New  York,  NY  10003 

14.  MONITORING  AGENCY  NAMF.  6 AODRESSfif  dllterent  tronx  Controlling  Office) 


13.  NUMBER  OP  PAGE^, 


15.  SECURITY  CLASS?75TfJ 


Unclassified 


I5a.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (o (title  Report) 


This  document  has  been  approved  for  public  release  and  sale 
distribution  is  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  l/io  abafrocf  entered  In  Blech  30,  It  different  from  Report) 


Rsprinted  from  AGARD  Conference  Proceedings  No. 
Subject:  Unsteady  Phenomenon  in  Turbo  Machinery 


18.  SUPPLEMENTARY  NOTES  I 


19.  KEY  WORDS  (Continue  on  reverae  aid a II  necaaaary  *nd  Identify  by  block  number) 

Aeroelasticity  Oscillating 

Nonstationary  Aerodynamics  Periodic  Flo 
Forced  Vibration  Flow  with  Ro 

Computation  Fluid  Dynamics 


AB^RACT  (Continue  on  reverie  eld*  //  neceoemry  en d identity  by  block  number) 


t^The  unsteady  response  of  a flat  plate  airfoil  to  large  non 
potential  flow  disturbances  is  evaluated  in  the  form  of  a 
translating  rectangular  grid  of  eddy-array.  The  singularity 
distribution  principle  in  combination  with  the  time-marching 
technique  was  used  as  the  method  of  analysis.  — 


FORM 
1 JAN  73 


EDITION  OF  I NOV  6S  IS  OBSOLETE 
S/N  010  7”  014*460  I | 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  Data  Rnlarad) 


